Magnetic field asymmetry of non-linear transport in Aharonov Bohm rings 
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Fundamental Casimir-Onsager symmetry rules for linear response do not apply to non linear trans- 
port. This motivates the investigation of nonlinear dc conductance of mesoscopic GaAs/GaAlAs 
rings in a 2 wire configuration. The second order current response to a potential bias is of partic- 
ular interest. It is related to the sensitivity of conductance fluctuations to this bias and contains 
information on electron interactions not included in the linear response. In contrast with the linear 
response which is a symmetric function of magnetic field we find that this second order response 
exhibits a field dependence which contains an antisymmetric part. We analyse the flux periodic and 
aperiodic components of this asymmetry and find that they only depend on the conductance of the 
rings which is varied by more than an order of magnitude. These results are in good agreement 
with recent theoretical predictions relating this asymmetric response to the electron interactions. 



X 



Magneto-transport properties of mesoscopic systems 
are sensitive to interferences between electronic waves 
and exhibit characteristic signatures of phase coherence. 
!l| Among them are universal conductance fluctuations 
(UCF) leading to reproducible sample specific magne- 
toresistance patterns, which in a ring geometry are mod- 
ulated by the flux periodic Aharonov-Bohm (AB) oscil- 
lations. Beside these effects on the linear conductance, 
it has been shown that mesoscopic systems exhibit rec- 
tifying properties related to the absence of spatial in- 
version symmetry of the disorder or confining potential. 
This gives rise to a quadratic term in the /, V relation: 
I = GiV + G2V 2 . This non linearity was predicted the- 
oretically [3] and observed experimentally more than 15 
years ago [3j. It was understood as a direct consequence 
of the sensitivity of conductance fluctuations to the Fermi 
energy with a characteristic scale given by the Thouless 
energy, E c — h/ro, where Trj = L 2 / D is the diffusion 
time across the sample of size L. More recently it has 
been pointed out that non linear transport coefficients do 
not obey Casimir Onsager symmetry rules Q and may 
include, in a 2 wires measurement, a component antisym- 
metric in magnetic field with a linear dependence at low 
field. On the experimental side the existence of a term 
linear both in field and current was found in macroscopic 
helical structures and attributed to magnetic self induc- 
tance effects [f| . It was subsequently observed in carbon 
nanotubes @ and suggested to be related to their helical 
structure. More generally, this field asymmetry of non 
linear transport has been theoretically shown to be re- 
lated to electron-electron interactions both in chaotic Q 
and diffusive [|[ systems. At the single impurity level, 
this effect can be simply viewed as the modification of 
the electron density dn(r) around a scatterer in the pres- 
ence of a current through the sample [§, Due to 



Coulomb interactions this results, in a modification of 
the potential around the impurity by a component lin- 
ear in current. In a phase coherent sample, this bias 
induced change of disorder potential dUdis modifies the 
conductance fluctuations giving rise to the nonlinear con- 
ductance G2 defined as: G(V) — G\ + G2V + ... Just 
like the chemical potential measured in a multiprobe 
transport measurement [l[, dUdis exhibits field depen- 
dent fluctuations which have both symmetric and anti- 
symmetric parts, including a component linear in B at 
low field. As a result G2 has a symmetric component 
in B and an antisymmetric one, respectively equal to 
G%' AS = (G 2 (B) ± G 2 (-B))/2 which both vary on the 
flux scale $ c characteristic of conductance fluctuations. 
The antisymmetric component only exists in the presence 
of electron-electron interactions. The typical amplitudes 
of these components have been calculated in a 2D system 
0)@)[IH and yield for weak interactions: 



6G% ~5Gi{e/E c ), SG^ S 



fint — f(\mc\)SGl (1) 

9 



where the function f(x) is equal to x for < x <C 1 and 
to 1 for x >■ 1. g —< G\ > and 5G\ ~ 1 are the average 
conductance and the typical amplitude of G\ fluctuations 
in units of e 2 /h. In the weak interaction regime investi- 
gated in [8] ji n t = 2vdUdis{'i f ) / dn{r) <C 1, where v is the 
density of states per unit surface. In the self-consistent 
treatment of Coulomb interaction [ll|, [ll|, the interac- 
tion constant is determined by the ratio of typical charg- 
ing energy of the sample ~ e 2 /C, and mean level spacing 
A as j int = 1/(1 + CA/2e 2 ) (the limit j int « 1 of Eq. 
□ corresponds to dll/dn = e 2 /2C)[l]|. The l/g factor in 
EqQ] indicates that the field asymmetry should be easily 
detectable in low conductance samples where g does not 
exceed 10, but is not observable in metallic mesoscopic 
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FIG. 1: Field dependence of Gf' AS and G%' AS ,(G AS shifted). 
These data were obtained on ringi in its initial state for a 
current of lOnA. Inset: micrograph of ringi. 

systems [3] ■ Carbon nanotubes [HI are in principle good 
candidates, but it is delicate to distinguish the effects due 
to the tube helicity from mesoscopic ones. Semiconduct- 
ing samples arc well suited for such investigations since 
they combine rather low conductance and large sensi- 
tivity to small fluxes [lj| [l?], [Hj]- We have measured 
the non linear conductance of 2 terminal GaAs/GaAlAs 
rings with only a few conducting channels. We find that 
Gi , like G\ , exhibits both AB oscillations and UCF con- 
ductance fluctuations. We show evidence of the existence 
of an asymmetry in magnetic field in G*2 from which we 
deduce the amplitude of the electron interactions. 

The square rings investigated in this experiment were 
obtained by shallow etching through an aluminum mask, 
of a 2D electron gas (2DEG) of density n e = 3.8 x 
10 15 m -2 at the interface of a GaAs/GaAlAs heterojunc- 
tion with Si donors. Contrary to previous experiments 
[3 [MM 03 > there is no electrostatic gate on the sam- 
ples [19[ • Due to depletion effects the real width of the 
rings is smaller than the etched one and is determined 
from magneto-transport data. We present data on 2 
rings, of circumference L — 4.8/ito and respective width 
W = 0.3, 0.45 ± 0.05/im. The elastic mean free path l e 
extracted from the conductance at 4K varies between 1 
and 2/ito, which is less than the value of the initial 2DEG 
and comparable with the side of the square ring. In situ 
modifications of the samples were obtained by short il- 
luminations with an electro-luminescent diode resulting 
in an increase of width and conductance of the rings. It 
was also possible to change the disorder configuration by 
applying current pulses in the 10 to 50 fiA range which 
decrease the average conductance of the samples. We 
could therefore, with only 2 samples, investigate a con- 
ductance ranging from g = 1 to 20. The measurements 
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FIG. 2: (a)Amplitudes of the UCF and AB components of 
Gi in units of e 2 /h, Gf and G AS in (fiV)" 1 as a function of 
the conductance of the rings. The open and closed symbols 
correspond to UCF and AB integrated peaks (straight lines 
are guides to the eyes). (b)Fourier transforms of G2 and Gi 
renormalised so that the amplitude of the AB peaks are iden- 
tical, (c) Conductance dependence of the ratio 5G AS /<5Gf for 
UCF (open stars) in comparison with Eq[2] (diamonds), where 
Jint = 0.85 for ringi and 0.94 for ring2. 



were conducted via filtered lines in a dilution refrigera- 
tor between 25mK and 1.2K. The samples were biased 
with ac current of frequency u>/2tt ~ 30Hz in the few nA 
range and voltage was measured with a low noise am- 
plifier followed by lock-in amplifiers detecting the first 
and second harmonics response V\ cos Lot and V2 cos2wi. 
The amplitude of the ac modulation was chosen to max- 
imize the second harmonics signal in the regime where it 
is still quadratic with the current modulation amplitude. 
In this regime the second order conductance G2 is simply 
related to V2 and V\ by G2 — ~ 2V2I /V^ . As shown in 
Figfflboth Gi{B) and G 2 {B) exhibit h/e periodic AB 
oscillations modulated by UCF fluctuations. 

More remarkable, whereas G±(B) is a symmetric func- 
tion of magnetic field as expected in a 2 wires config- 
uration, G2(B) exhibits a component antisymmetric in 
field G^ 3 ■ To compare the field asymmetry of G2 on 
the various samples (Fig. 2), we extract the amplitude 
of UCF and AB oscillations from the Fourier transform 
of G^ AS {B) and G X {B). The integrated UCF and AB 
peaks are noted 5G2' AS and <5Gf' AS and similarly for 
G\. These averages performed on a flux range much 
larger than <!> c do not depend on this range. We find 
that all these quantities are aligned on logarithmic plots 
as a function of g. We first note the relatively large ampli- 
tudes of the AB oscillations dG^ ' AS which are of the order 
of the UCF components SG^ ' S in contrast with the re- 
lated quantities in G\. The decrease of 8G\ and 5G\ to 
values much below 1 at low g seems to be at odd with 
the universal character of conductance fluctuations in G\ 
established deep in the diffusive regime j > 1. However 
this universal regime is not expected for g ~ 1. The even 




FIG. 3: (a)Temperature dependence of the amplitude of the 
AB oscillations in Gi, Gf and (ring2). Continuous lines 
are the fits corresponding to the regions (1) &bT <C E c (2) 
fcsT ^> E c but L <C L<f, (3) L > L^>.(b)Examples of tempera- 
ture dependence of G2 in ring2 for different values of magnetic 
held between -1000 and 1000 Gauss. 

larger dependence measured in <5Gf (<?) indicates, accord- 
ing to ([1]), that E c — gA decreases with g. This finding 
can be attributed to the physics involved in the illumina- 
tion process which increases the width of the rings (re- 
sulting in a decrease of A) and also decreases the elastic 
scattering time. Finally we find that the ratio 5G$ S /5G% 
depends only slightly on g and is equal to 0.3 ± 0.2(see 
Fig. 2c). This result is a priori in contradiction with EqQ] 
from which a 1/g dependence of 8G% /8G§ is expected. 
However Eq[T]is only valid for 7^ <C 1 in partially closed 
quantum dots, i.e. whose classical resistance is domi- 
nated by the contacts but still much lower than quantum 
resistance. It was recently found by one of us [20| using 
random matrix theory that in these systems (5Gf also 
strongly depends on 7, n t, when it is not negligible com- 
pared to 1, and in the limit $ 3> $ c : 

5Gi s /SGl = l/ v /l + 2(.9/JG 1 ) 2 (l/ 7m t-l) 2 (2) 

identical to ([T]) in the limit of small 7^ . It is interesting 
to note that EqfJ] predicts <5G^ S /<5Gf independent of g 
in the limit ji nt = 1. It describes remarkably well our 
experimental results, yielding 7 i nt ~ 0.94 ±0.02 for ring2 
and 7 i„t ~ 0.85 ± 0.02 for ringi . These values are very 
close to the estimated values of 7 j n j = 1/(1 + GA/2e 2 ) ~ 
0.98 from the geometry of our samples. Our determina- 
tion of the interaction constant from the dimension less 
quantity 5Gf s /SG^ is much more accurate than the es- 
timation done in reference [17j | from the analysis of the 
low field component of Gf s only, but consistent with it. 
Finally, our results show that our rings are in a highly in- 
teracting regime and behave as partially closed quantum 
dots due to their relatively long mean free path. 

~ ~ S AS 

The temperature dependence of 5G\ and 5G 2 ' are 
shown on Fig(3^,. As already observed [22j 5G\ is only 
weakly T dependent below the Thouless energy like 
1 - AT 2 with A = 2(k B /E c ) 2 and decays at higher tem- 
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FIG. 4: (a)Low field dependence of G2 on ringi for differ- 
ent temperatures, (b) Field dependence of G;f S after low pass 
(light line) and high-pass (bold line) filtering , on ringi and 
ring2 respectively lower and upper curves (shifted for clarity). 
Note the extinction of AB oscillations in the vicinity of B — 0. 



perature as T Deviations above 0.5K are consis- 

tent with exp(— L/L^T)) with a T dependence of the 
phase coherence length as T" 1 / 2 but this last fit on a 
small range of temperature is not unique and just in- 
dicative. SG2 ' AS have similar T dependences, which are 
nearly identical to (SGi(T)) 2 with a T" 1 decay in the 
limit ksT 3> E c with L C in agreement with theo- 
retical predictions [H, [ll[ . The temperature dependence 
of G2 at fixed magnetic fields is depicted in Fig[3)D. In the 
same way as observed on Gi [23[, Gi(T) exhibits a non 
monotonous variation with temperature on the scale E c 
which randomly fluctuates with magnetic field. In some 
cases we also observed (see FigH^,) that the phase of the 
AB oscillations in G2 depends on temperature. Surpris- 
ingly it remains pinned either to or 7r at zero field with 
the appearance of a second harmonics contribution in the 
region of temperature where the sign change occurs. At 
larger field the phase takes any value between and it. 
Note that this last effect is observed both in Gi and G2 
as a result of phase modulation of the AB oscillations 
by the UCF, but these phase modulations are symmetric 
in field on Gi, and not on G2. In short, we find that 
AB oscillations on G2 are symmetric in zero field and 
the asymmetry only appears at higher field. This is also 
clearly seen in the AB oscillations on Gf s . We observe 
in all samples that their amplitude vanish linearly at zero 
field (Fig.4b). 

We now propose a simple explanation of the extinc- 
tion of AB oscillations at low field in G^ S linked to 
the larger AB component observed in G2 compared to 
G\. In the semiclassical approximation (kfl e 3> 1 ) it is 
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possible to express the conductance of a phase coherent 
ring in terms of interference between scattering ampli- 
tudes of electronic waves, G = Re [EijA^Aj expi4>ij(B)], 
where indices i and j run on all pairs of possible dif- 
fusive trajectories going from one terminal to the other 
and 4>ij(B) = 4>i — <fij + 2irBSij /$o- The phase at zero 
field fa is the integral (1/h) f. Udi S {r(t))dt on the diffusive 
trajectory i assuming that the screened disordered po- 
tential Udis varies smoothly on the electron Fermi wave- 
length which is a reasonable assumption in a 2DEG. Sij 
is the surface comprised between trajectories i and j. On- 
sager symmetry rules imply that the 2 terminal magneto- 
conductance takes the following form at zero tempera- 
ture: 

G Y.,.,A,A, cos(<f>ij)cos(2irBSij/$ ) (3) 

In the presence of a current through the sample asso- 
ciated to a potential drop V, the local electronic den- 
sity and consequently the scattering potential is modi- 
fied with a term dUdis(r,V) which induces phase shifts, 

d &j = Ji~ fjil/tydUdisiritydt. The quantity G 2 = 
(dG(V) j dV)v=o is then directly related to these phase 
shifts through: 

G 2 = XijAiAjidfaj/dV) sin(^) cos(2nB S i:j /$ ) (4) 

Since d(j)ij increases with the length of the interfering tra- 
jectories i and j, long trajectories encircling the ring (AB 
oscillations) contribute more than trajectories within the 
same branch of the ring (UCF). This explains the larger 
relative amplitude of the AB oscillations and the larger 
harmonics content in G 2 compared to G\ (see Fig. 2 a,b). 
In a diffusive open system dUdis is of the order of eV 
and the main contribution to the typical value of d4>ij is 
eVTo/h and is independent of B leads to the expression 
of <5G 2 of EqfT] This contribution is strongly attenuated 
for a ballistic system with strong interactions, like in [llj ]. 
where most of the potential drop takes place at the con- 
tacts. Other contributions to <5G 2 of the order of 1/g 
are obtained by taking into account mesoscopic non local 
field dependent fluctuations in the potential dUdi S (r, B), 
with symmetric and antisymmetric components in B, 



both of order of jmteVSGi/g at high flux compared to 
$ c . The antisymmetric component in dUdi S (r, B) is ex- 
pected to typically vary like Jint(§Gi/g)V<&/$ c at low 
flux. We can then deduce from Eqf5]the main contribu- 
tions to the AB oscillating part of G^ 5 , which originates 
from pairs of trajectories encircling the rings: 

G$ s = „.,>:,..,< .l,.l ? r, , sin(4 3 )^- cos(2tt$/$ ) 

(5) 

There is also a smaller contribution of terms in 
sin(27r$/$o) of the order of SGi/SGi which is less than 
0.2 in our experiments. This provides an explanation for 
the linear increase of the amplitude of the AB oscillations 
of g£ s in a flux range <i> c for rings with g > 1 as shown 
in Fig. 4b. However this effect may not exist in very nar- 
row rings [HJ where conductance fluctuations are not ob- 
served. It is interesting that this simple heuristic model 
can also explain the larger value of G 2 in a diffusive com- 
pared to a ballistic system and the different values of the 
ratio G^/Gf of the order of Jint/9, respectively 1, in a 
diffusive, respectively ballistic, system. 

In conclusion we have shown evidence of a field asym- 
metry on the second order response of GaAs/GaAlAs 
rings of mesoscopic origin which contains both AB os- 
cillations and conductance fluctuations. This asymme- 
try is characterized by JG^/JGf and analyzed within 
theoretical predictions expressing this ratio with only 2 
parameters, the dimensionless conductance of the rings 
and the interaction constant whose value can be deter- 
mined ji n t — 0.90 ± 0.05. We have also found that the 
relative amplitude of the AB oscillations compared to the 
UCF is much larger in G 2 than in G\ with the existence 
of a linear low field modulation in the AB oscillations in 
the antisymmetric component of G 2 . These effects can 
be understood within a simple semi-classical description 
of quantum interference. We thank D. Mailly and F. 
Pierre for their helpful contribution to the fabrication of 
the samples as well as M. Biittiker, M.Ferrier, G. Mon- 
tambaux, B.Reulet, B. Spivak and C. Texier for fruitful 
discussions. Note added During completion of this work, 
related experimental work on gated quantum dots (l7j 
and small Aharonov Bohm rings [1 81 ] appeared. 
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